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Abstract — Consider the problem of joint uplink scheduling 
and power allocation. Being inherent to almost any wireless 
system, this resource allocation problem has received extensive 
attention. Yet, most common techniques either adopt classical 
power control, in which users are received with the same 
Signal to Interference-plus-Noise Ratio (SINR), or use centralized 
schemes, in which base stations coordinate their allocations. 
Modern systems such as IEEE 802.16e/m or 3GPP LTE, however, 
require inter-cell interference mitigation without compromising 
the benefits of distributed algorithms. 

We suggest a novel scheduling approach in which each base 
station, besides allocating the time and frequency resources 
according to given constraints, also manages its uplink power 
budget such that the aggregate interference, "Noise Rise", caused 
by its subscribers at the neighboring cells is bounded. Our 
suggested scheme is completely distributed between the cells, 
requiring neither coordination nor message exchange. 

We rigorously define this allocation problem under noise rise 
constraints, give the optimal solution and derive an efficient 
iterative algorithm for it. We then discus a relaxed problem, 
where the noise rise is constraint separately for each sub- 
channel or resource unit. While sub-optimal, this view renders 
the scheduling and power allocation problems separate, yielding 
an even simpler and more efficient solution, while the essence of 
the scheme is kept. Via extensive simulations, we show that the 
suggested approach increases overall throughput dramatically, 
with the same level of fairness and power consumption. 

I. Introduction 

The desire to provide integrated broadband services while 
maintaining Quality of Service (QoS) guarantees bestows 
growing interest in schedule access techniques used in 
multiple-access protocols for future broadband radio sys- 
tems. Hence, for example, IEEE S02.16e,m (WiMAX) (2, 
[|2] adopted Orthogonal Frequency-Division Multiple Access 
(OFDMA) as the access mechanism. The design of radio 
resource mechanisms, and in particular scheduling and power 
control, while maintaining some type of fairness, is challeng- 
ing due to the high variability in link-quality. In this work, 
we consider the joint uplink scheduling and power control 
problem for wireless networks. 

Common power-control approaches to the uplink resource 
allocation problem are either to assign transmission power to 
the mobile stations such that all users are received at the base 
station with the same Signal to Interference-plus-Noise Ratio 
(SINR) 0, 01, or to allow users to transmit at their maximal 
available power [5|. Both of these techniques optimize the 



user/cell throughput, neglecting the interference injected to 
neighboring cells. However, since OFDMA systems, for ex- 
ample, are sensitive to inter-cell interference, the interference 
from neighboring cells can dramatically decrease the SINR 
received at the base station, hence reduce the mobile station 
throughput. As a result, the role of the power control becomes 
decisive in providing the required SINR while controlling at 
the same time the interference caused to neighboring cells. 
This is challenging, as on the one hand a mobile station near 
the base station is expected to have high quality link, hence 
high throughput, even when transmitting in low power, while 
a distant mobile needs to transmit in a much higher power 
to attain the same throughput, and on the other hand as far 
as inter-cell interference is concerned, mobiles near a base 
station can transmit in high power since they are not in the 
proximity of other cells, while distant mobiles which can be in 
the proximity of other cells should not transmit in high power 
as they can interfere with other (neighboring) base stations. 

In order to limit the interference to neighboring cells, 3GPP 
has approved the use of Fractional Power Control (FPC) [6|. 
This new proposal makes users with a higher path-loss operate 
at a lower SINR requirement so that they will more likely gen- 
erate less interference to neighboring cells. According to this 
approach, to preserve fairness, most of the resources should 
be allocated to far-away users. Similarly, the IEEE 802.16m 
power control scheme deducts a fraction of the downlink 
signal-to-interference ratio (SIR) from the transmission power. 
By that, it reduces the interference caused by cell-edge users. 
The 3GPP LTE standard also suggests a different approach for 
inter-cell interference handling termed Inter-cell interference 
coordination (ICIC) (see, e.g., J7J). ICIC provides tools for dy- 
namic inter-cell-interference coordination of the scheduling in 
neighboring cells such that cell-edge users in different cells are 
preferably scheduled in complementary parts of the spectrum 
when required. However, ICIC requires coordination between 
neighboring cells both in terms of exchanging information 
regarding subscribers at one cell and their interference level on 
other neighboring cells, as well as coordination in the resource 
allocation, which further complicates the scheduling process. 

In this work, our contribution is three-fold. First, we in- 
troduce a different approach, which controls the inter-cell 
interference, yet does not require any cross-deployment com- 
munication or coordination. In our approach, the uplink inter- 



2 



cell interference that each base station is allowed to create 
is bounded. This limited egress interference budget, termed 
Noise Rise, is then allocated to users in conjunction with the 
ordinary resources (time and frequency) allocation, according 
to some fairness criterion and channel condition. Controlling 
the average interference level sensed by each base station 
provides a more predictable uplink SINR, which allows a 
lower interference margins and a more efficient rate selection. 
Hence, it obtains a higher capacity and a better coverage. 

Second, we formalize the scheduling problem under the 
noise rise constraint as a constraint convex optimization prob- 
lem, and provide an efficient iterative algorithm that is proved 
to solve it optimally. Moreover, we suggest a second setting, 
in which instead of bounding the average noise rise over all 
the channels, allowing some sub-channels to contribute more 
noise rise at the expense of further limiting the noise rise on 
others, we bound the noise rise on each sub-channel to the 
exact same value. The later settings allows the decoupling 
of the scheduling algorithm from the power control and thus 
facilitates an even simpler algorithm. 

Third, we thoroughly evaluate the noise rise concept via an 
extensive set of simulations, using both the direct analytical 
setting defined by the optimization problem and a more 
elaborate, all-inclusive setting defined by IMT- Advanced |8|. 
These simulations clearly depict that the suggested approach 
dramatically increases the overall throughput achieved in each 
cell compared to the traditional approach, while maintaining 
a high level of fairness. 

II. Noise Rise 

Traditional cellular communications such as CDMA net- 
works suffer from intra-cell interference caused by the pseudo 
orthogonality of the intra cell CDMA codes. Alternatively, the 
orthogonality of the sub-carriers in OFDMA networks solves 
the intra-cell interference resulting in inter-cell interference 
limitations. In such networks, the interference levels seen in a 
base station depends on the power control and scheduling of 
neighboring cells and not only on its own resource allocation. 
Accordingly, an optimal scheduling and power allocation 
would require coordination and synchronization across all 
deployment. Clearly, such a synchronization could be pro- 
hibitively complex and impractical. Without cross deployment 
coordination, where each cells performs its power control and 
scheduling independently of its neighbors, the interference 
profile seen in uplink transmission becomes highly dynamic. 
To overcome the variability of the inter-cell interference a 
high interference margin is accounted for while selecting the 
transmission rate and coding scheme. That is, transmitting in 
a more robust rate (lower rate) to maintain reliability at worse 
cases of the interference pattern. Clearly, a hight noise margin 
assumption reduces the cells uplink spectral efficiency. 

In addition to the cell capacity, constrained interference 
is vital for cell coverage. If the noise plus interference is 
not constrained, due to the maximal power transmission of a 
mobile station, certain users (typically at the cell edge) would 
not be able to maintain a reasonable communication capacity. 



These users would be blocked and removed from the coverage 
area, resulting in a reduced cell coverage. 

Typically, the interference is represented as a measure of 
noise rise, which is defined as follows. 

Definition 1. Uplink noise rise is defined as the total uplink 
received noise plus interference power over the background 
noise power. Formally, let 7 denote the noise rise, then Nq + 
7 = 7- No, where I is the total interference at the base station 
(receiver), and Nq is the background noise. 

The constraint noise rise approach, which is described in the 
next section, provides a fully distributed mechanism (without 
a cross deployment coordination and synchronization) that 
dramatically reduce the variability of the noise rise (due 
to inter-cell interference) allowing a more aggressive rate 
selection (taking a much smaller noise rise margin). 



A. Constrained noise rise approach 



The constrained noise rise approach aims at allowing the 
scheduler to operate in a fixed Noise plus Interference condi- 
tions and facilitating a more aggressive rate selection. Specif- 
ically, we aim at bounding the egress interference and show 
that it also bounds the noise rise seen by uplink transmission 
at all cells. Formally, let ftA{k) denote the set of backlogged 
users at base station k. Denote by Q(k*) the set of neighboring 
base stations of base station k* and by Ii n (k*) the ingress 
interference level per resource unit at base station k* receiver, 

hn{k*)= ^ L i(k),k' ■ Pi(k) (1) 

keG(k') i(k)eM(k) 

where, L^k).k' denotes the channel gain between user i(k) 
of base station k and base station k* , and Pi(k) the transmission 
power of user i(k). In fl}, the first summation is over all 
neighboring cells except the serving cell k*, and the second 
summation is over all users in the neighboring cells. 

On the other hand, the egress interference by the users of 
cell k* to the neighbors is given by 

Ieg{k*) = ^ ^ Li(k').k ■ Pi(k*)- (2) 

i(k*)<£M{k') fceS(fc*) 

Assuming a fully homogeneous and symmetric deployment, 
we have that on the average the ingress interference that 
neighboring cells injects to a base station equals to the egress 
interference caused by the base station to the neighboring cells. 
That is, 
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(3) 



Therefore, maintaining a constrained egress interference at 
all deployment, realizes our goal of a constrained ingress 
interference. 

Let li denote the normalized interference that a user i's 
uplink transmission injects to its neighboring cells, that is 

k = E £i(fc*),fc- (4) 
keg(k*) 

The interference caused by the users of cell k*,I eg (k*) , can 
be written as 



Ieg(k*) — 2J k(k») ■ Pi(k*) 

i(k*)eM(k') 



(5) 



We are now ready to write the noise rise constraint on the 
scheduler and power control allocation. The constraint bounds 
the egress interference caused by the cell uplink transmission 
such that the ingress interference at all the deployment would 
remain bounded as well. The scheduler and power control 
should comply with the following constraint: 



i(k')GM(k*) 



h(k*) ' Pi{k') < I 



(6) 



We note that (O implies that the scheduler may consider 
noise rise as a resource to be allocated to transmitting users. 
That is, instead of considering the traditional power allocation, 
the noise rise approach considers noise rise allocation, i.e., 
li -pi. Accordingly, for example, if the egress interference of a 
given user is high, the scheduler should reduce its transmission 
power and increase its allocated bandwidth. Alternatively, the 
schedule could co-allocate this user with another user that 
"consumes" less noise-rise such that the noise rise budget is 
kept. 

In the following sub-section we provide the means for a 
base station to estimate the required normalized interference. 



B. Normalized interference 

To comply with (O each base station should have informa- 
tion on the normalized interference Typically, a base station 
cannot directly measure this coefficient. In fact, the normalized 
interference can be measured at the surrounding cells and 
reported to the base station (e.g., via Inter Cell Interference 
Coordination). Alternatively, a base station may estimate the 
normalized interference from its users' downlink channel 
state reports without the need for inter-cell communication 
or coordination. 

The downlink SIR (signal-to-interference ratio) measured 
by a user is given by 



SIR?(k' ) 



T PDL 

^k*,i(k")" 

E L kAk , } PDL 

fces(fc*) 



(7) 



where P DL is the base station downlink transmission power. 
With the channel reciprocity [9|, i.e., L^.i — L; L , k , we have 

Li, 



H(k*),k" 



keG(k*) 



i(k*),k 



Therefore, 



fceS(fc') 



^i(k*),k* 



(8) 



(9) 



Finally, we note that in an interference-limited scenario (where 
the interference is much larger than the background noise), 
which is typical, one could use the downlink SINR instead of 
the downlink SIR, neglecting the background noise. 

III. System Model 

We consider the problem of scheduling and resource al- 
location for the uplink of an OFDMA cell, where a set 
M. = {1, ...,M} of users (with backlogged traffic) transmit to 
the same base station. The total frequency band is divided into 
a set TV = {1, N} of distributed resource allocation units 
(e.g., frequency sub-bands or sub-channels). The distributed 
resource allocation units are of the same size as a physical 
resource units that has undergone sub-band partitioning and 
permutations. The permutation spreads the resource units' sub- 
carriers across the whole frequency band. Accordingly, for 
each user i the channel quality of all logical resource units 
can be assumed identical. The channel quality of a resource 
unit represents a collective quality indicator for the resource 
unit, e.g., the average across the sub-carriers in the unit. 

Let Xi, < Xi < 1, be the fraction of the resource (i.e., 
fraction of the frequency band) allocated to user i. The total 
allocation across all users should be no larger than 1, i.e. 



E 



Xi < 1 



(10) 



In practical systems, only a single user can transmit over a 
given resource unit (assuming a single user MIMO), accord- 
ingly, Xi is a fraction of the total number of resource units. We 
consider systems with a large number of resource units where 
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Xi could take any fractional value neglecting the rounding 
errors (e.g., 92 resource units for a 20MHz band in IEEE 
802.16m). Note that the rounding here differs from that of [5], 
which considers the fraction of a single resource unit (i.e., Xi 
could take only or 1). Differently, we consider the fraction 
of the divisible whole frequency band. 

Denote by pi the transmission power of a user i over 
all resource units, which is subject to a maximum power 
level, pi < P. In this paper we consider interference limited 
scenarios where the transmission power is bounded by the 
constrained noise rise rather than the mobile transmission 
power capabilities. Accordingly, in the following we shall 
assume that users do not have a maximal power constraint. 

Time is divided into equal length time slots (i.e., sub-frames 
in the IEEE 802.16m terminology). At the beginning of every 
time slot, the scheduler seeks to maximize a (time-varying) 
weighted sum of the users rates. 

The scheduler and power control mechanism are assumed to 
have knowledge of the channel gain, which comprises of the 
long term parameters of the link between the base station and 
the user, such as path loss and shadowing factor, as well as the 
short-term time-varying spatial fading. Let Li be the channel 
gain from user i to the serving base station. Additionally, the 
scheduler is aware of the target Interference plus Noise at the 
serving base station, which translates to a normalized received 
signal-to-interference plus noise ratio (SINR) per unit transmit 
power, ej , for each user i. 

We assume that n (t) is the achievable rate of user i at time 
t in a Gaussian Multiple Access Channel using timesharing 
[ 1 1 . Formally, ri(t) is given by 



n(t) = Xilog l + 



where the tuple (x,p) S x denotes the resource and power 
allocation and the set \ is given by 

X = \ (x,p) > : x i < !. Yl **(*) P^ 1 ] (12) 

I ieM ieM ) 

k(t) being the normalized interference of user i at time 
slot t. We use bold symbols to denote vectors, e.g., x = 
{xi,x 2 ,...,x M }, P = {pi,P2,---,Pm}- We adopt the 
gradient-based scheduling framework lfTTl - lfT5l . Specifically, 
the scheduler solves the following optimization 



max uji(t)ri(t) 



(13) 



ieM 



where uii (i) > is a time-varying QoS weight assigned to the 
i'th user at time t. 

We concentrate on weights that depend on the average 
throughput attained by each user up to the t-th slot, and capture 
some fairness notation. For example, u>i(t) = jrjjji where 
Ti(t) is the average throughput of user i at time t, which 
captures proportional fairness ifTTl . Ifl4l . fl5l . 

Note that ( fT3l must be re-solved at each scheduling instant 
(e.g., each sub-frame) because of changes in both the resource 



unit state as well as the weights. 

For the ease of presentation, in the following we shall omit 
the time index t. 

IV. Joint uplink scheduling and power control 

At the beginning of each time slot, the base station schedules 
a subset of the backlogged users to the available resource units 
and assigns transmission power to each scheduled user. The 
base station aims at maximizing the achievable rate while 
providing users with a fair share of resources (according 
to a predefined fairness metric) and maintaining a bounded 
interference with neighboring cells. 

In conjunction with the scheduler and the power con- 
trol, a rate adaptation mechanism adjusts the transmission 
rate according to the allocated power. Clearly, the allocated 
throughput in a time slot (sub-frame) to a scheduled user is 
derived from the allocated resource units and the allocated 
rate (i.e., power). In other words, the scheduling and power 
allocation are coupled, in devising the allocated throughput, 
and should be performed jointly. The power control scheme 
optimizes the trade-off between allocated rate and contribution 
to the overall Noise Rise. Typically, users that are far from the 
base station are required to transmit in a high power in order 
to maintain a reasonable rate. However, these users are closer 
to neighboring cells, hence contribute more to the noise rise 
in those cells. That is, on the other hand, cell edge users are 
required to transmit in a low power in order to bound the 
interference with neighboring cells. 

At the beginning of each time slot the scheduling and power 
control scheme selects a feasible resource and power allocation 
tuple (x,p) that complies with the noise rise constraints and 
maximizes a time-varying weight assigned to each user, i.e., 



maximize y. u i ' r i(xi,Pi) 
(x ' p)e * ieM 



(14) 



where Ti{xi,pi) is the rate related to the resource and power 
allocation, and uii > is the time-varying weight assigned to 
the i-th user at the beginning of the time slot. These weights 
are the gradient of an increasing concave utility function of 
each user. Taking the rate as ri(xi,pi) — Xi log ^1 -I- ^r^j 
(|TTJ, we can formulate the joint power control and scheduling 
with noise rise constraint optimization problem. 

For each time slot, find the channel allocated to each 
user, denoted by by a; = {x±, x%, . . . , Xm}, as well as 
the set of power assigned to each user, denoted by p = 
{p 1 ,p 2 ,-.- ,Pm}, such that 



maximize 

sc,p 



1 °g( 1 + v 1 ) 

subject to Xi > 0, Vi G M, 

p l > 0, Vi S M, 



(15) 



l *Pi ^ L 



ieM 



5 



Proposition 1. Optimization problem H15j is convex with 
linear constraints. 

For a proof, see Appendix lAl 

From Proposition Q] it is clear that the optimal solution can 
be found numerically using standard optimization techniques. 
For example, it can be found through a similar method to 
that in J5), |[T6l . Yet, this direct approach might be pro- 
hibitively complex and infeasible for practical implementation 
in commercial base stations. Accordingly, a simpler solution 
is called for. In the next section, we give an efficient iterative 
algorithm, which uses the analytical solutions to two related 
sub-problems, to solve the above problem. Moreover, we show 
that the iterative algorithm converges to the global optimum. 



V. An Analytic Solution to the Joint Scheduling 
and Power Control Problem 

In this section we show that the optimal solution can be 
viewed as two intertwined water-filling-like problems, facil- 
itating a highly efficient solution which solves the complete 
optimization problem by fixing a subset of the variables each 
time (either powers or bandwidth) and solving the resulting 
watter-filling problem. We show that this iterative procedure 
is bound to converge, and, moreover, give analytical bounds 
on the possible values of the slack variables in each of the 
separate water filling problems, enabling us to converge to 
their solutions using a fast binary search. 

The optimization problem we discuss is as follows. 



maximize 

x,p 



ujiXi log | H 

Xi 



ieM 



subject to Xi > 0, Vt G M, 

Pi > 0, Vi e M, 
E Xi = 



(16) 



ieM 



ieM 



Denote maxja:, 0} by 



We first consider the analytical 



solution to this problem. Proposition [2] below gives a set of 
equations which the optimal bandwidth and power allocations 
satisfy. 

Proposition 2. Consider the joint power and bandwidth opti- 
mization problem in \16\ . The optimal power and bandwidth 
allocations, {p*}ieM< an d {x*}ieM> respectively, satisfy 



1 



1 + 



where Xi is the solution to the following equation 

p*a p*eiUi 



uji log(l 



+ X-2 







(17) 
(18) 

(19) 



and Ai, A2 are set such that 

2> 



(20) 
(21) 



For a proof, see Appendix [ 



A. An Iterative Algorithm 

While Proposition |2] gives necessary and sufficient condi- 
tions for power allocations p and bandwidth allocations x 
to be optimal, its direct computation is cumbersome, as the 
equations for both types of variables are intertwined. However, 
in this sub-section, we show that the problem in dT6l > can 
be solved optimally by a highly efficient iterative algorithm, 
which, unlike standard iterative optimization procedures, does 
not jointly optimize on all variables, but rather utilizes the fact 
that when separating the power variables from the bandwidth 
ones, each optimization problem has a relatively easy watter- 
filling-like analytical solution. 

The important observation is as follows. Fixing the band- 
width variables {xi}t e M^ the resulting optimization problem 
is 



maximize 
p 



ieM 



ujiXi log [ 1 + 



subject to pi > 0, Vz e M, 
^ kPi = I- 

ieM 



(22) 



The solution to this problem is the well-known water filling, 
e.g., 1 17 Example 5.2]. Hence, it is easily solvable (note 
that the weights ujiXi and the noise rise constraints U only 
serve as scaling factors, and do not change the essence of the 
problem in the separated case). Fixing the power variables 
{Pi]ieM results in a relatively similar optimization problem, 
which although involves an implicit equation on each Xi, is 
also straightforward to solve. The iterative algorithm will then 
alternate between the two solutions, fixing one set of variables 
based on the results of the previous iteration. A pseudo code 
of the algorithm follows. 

Iterative- Water-Filling^, I, I) 
1 x <- x° : such that x° > Vi, Yh x i = 1 
repeat 



2 
3 

4 
5 
6 
7 



Ai 



SOLVE^ kXi 

J_ 

Alii 



Vi 



I) 



X 2 ,x 4- SolveQJj Xi = 1 



; log(l + aa 



until Converge 



Xi~\-piC 



e-+A 2 =o) 



8 return {x,p} 

When evaluating Algorithm Iterative-Water-Filling, the 
two key aspects are complexity and convergence. First, con- 
sider the number of operations in each iteration. The first step, 
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as mentioned, is a basic water-filling procedure. The value 
of Ai can be calculated by first ordering the users according 
to , then including users from the smallest to the largest 
until the "water level" 4~ satisfies the noise rise constraint. 
This is done in 0(M log A/). As for the second step, the 
solution is more involved, since it cannot be solved explicitly. 
However, as the following proposition states, the solution is 
monotonic in A2, with upper and lower bounds on the value of 
the optimal A2, hence can be solved efficiently by a logarithmic 
time binary search. 

Proposition 3. For each i G M, let Xi be the solutions to 



uj 1 log 1 



Pi^i 



J 7 



Pi&i 



+ A 2 



0. 



(23) 



Then, for all i G A4, every A 2 < and any LUi > 0,Pi > 
and ei > 0, we have: 

1) Xi > and Y2ieM ^ * s mon °t° n i ca ^y increasing in A2. 

2) The value of A2 in ( 1231 ) such that Xi — 1 satisfies 



A' 2 " m < A 2 < Xf ax , where 



\ mtn 

A 2 — 



mm <, uji 

ieM 

\ max 
A 2 



1 + Mpiei 



max <, uji 

ieM 



Pi&i 



1 



log(l + Mp,e, 



log(l +piei) 



Thus, when solving the bandwidth iteration in Algo- 
rithm Iterative-Water-Filling, simply compute A™ in 
and \™ ax , and apply a binary search for the value of A 2 such 
that I J2i x i~ 1 — e - The computational cost is 0(Af log(l/e)) 
assuming d23l is solved for Xi in 0(1) for fixed A 2 . 

For a proof, see Appendix Icl 

Finally, we mention two important results on the conver- 
gence of Algorithm ITERATIVE- Water-Filling. As men- 
tioned, the algorithm iteratively solves two optimization prob- 
lems, each one involving half of the parameters to be opti- 
mized (either powers of bandwidth allocations). This is an 
alternating optimization procedure. While this procedure can 
fail for some utility function^ it is important to note that in 
the specific case of our joint power and bandwidth scheduling 
with noise-rise constraint, it is bound to converge. 

Corollary 1. Assume Algorithm ITERATIVE- Water-Filling 
is used to solve ( 1761 ). Then the power assignments p and 
channel allocations x converge to the global optimum of U6h 

For a proof, see Appendix iDl 

Note that the result in [18 Theorem 2] implies that as 
long as the power and bandwidth allocations p, x are not in 
the range of the fixed points of T[p,x] (an iteration of the 
algorithm), the negative utility function is strictly decreasing 
in each iteration. Moreover, the result in Corollary Q] can be 
made even stronger if the negative utility function is strictly 

'E.g., x 2 — 3xy + y 1 when alternating between the optimization on y and 
the optimization on x 



convex in the range. Consider, for example, the two user case. 
While the function 

-wixi log H uj 2 x 2 log H I 



Xl 



X2 J 



is not strictly convex for all p\,pi, x\,x% (the Hessian matrix 
is not of full rank and hence not positive definite), it is strictly 
convex under the constraints x\ + x 2 = 1 and l±pi + I2P2 = I 
(for positive Xi and p{). In this case, by lfl8l Theorem 3], 
the alternating optimization in Algorithm Iterative-Water- 
FlLLING converges q-linearly to the global optimum, from 
any starting point (pi,p 2 , x\, a; 2 )° in the range. That is, each 
iteration of the algorithm decreases the distance to the global 
optimum by a constant multiplicative factor q € [0, 1). In other 
words, 

||(p,x)t-(p*,x*)||<q||(p,x)t-i-(p*,x*)||. 
VI. Constrained noise rise density 

In this section we consider a second noise rise settings 
where the noise rise density is constrained. That is, instead 
of noise rise budget over the whole frequency band, the noise 
rise per resource unit is constrained. 

The constrained noise rise density approach is also suitable 
to the contiguous resource allocation scheme (i.e., without any 
sub-band partitioning and permutations). It is also viable to the 
LTE SC-FDMA uplink allocation. 

Formally, the constraints on the noise rise density for each 
allocated user i is given by 



U- — <I Vi 



(24) 



It is easy to see, by summing both sides of d24l > over all i G M, 
that a resource and power allocation {x,p} that complies with 
d24l) will not exceed the noise rise constraint over the whole 
frequency band. 

The scheduling and power control problem with constrained 
noise rise density is as follows. 



Pie, 



maximize \^ uii{t)Xi log ( 1 

x.n — ' \ 

ieM v 

subject to xi > 0, Vi e M, 

Pi > 0, Mi G M, (25) 
^2 ^ < 1, 

ieM 

li ■ — < I, Vi : Xi > 0. 

Xi 

It is easy to see that bounding the noise rise density allows 
us to decouple the scheduling from the power allocation. That 
is, the power allocation and (consequently) rate assignment 
are determined by the required noise rise density regardless 
of the allocated resource units. Accordingly, for each user 
i the allocated power density (e.g., power per sub-carrier) 
is given by I 1 = j-, and the achievable rate of user i 

is given by Xi log ( 1 + ^ j . Now, let the user i* be such 
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that i* = arg max |tJi(i) log ^1 4- i. The power and 
frequency allocation that solves (l25l l. allocates all frequency 
band to user i* with the power p,* = J-. A pseudo code that 
depicts the algorithm follows. 

Resource- Allocation^, L, SIR dl ) 

1 X <r- 0,p <r~ 

2 for i e M 

3 do uji *r- Fairness -Weight(^) 

4 i. j L± 

5 

Xi Li 

6 n <- Rate-Adaptation(^) 

7 i* <r- arg max {wi • r A 

8 «— 1 

10 return {a;,p} 

Z, = {ii, L2, • • • j Lm}, where Li is the channel gain 
between user i and its serving base station. SIR DL = 
{SIR? L ,SIR% L , SIR^ L }, where SIR? L is the down- 
link signal-to-interference ratio at the receiver of user i. k 
is the user i normalized interference, uji is obtained through 
Fairness-Weight(£), which computes the user i fairness 
weight at time t. Typically, = jq-fj, where Ti(t) is the 
average throughput of user i at time t. 

Note that, the algorithm RESOURCE-ALLOCATION does 
not necessarily assume a Gaussian Multiple Access Channel 
(i.e., Ti — log ^1 + Instead, the algorithm calls a rate 

adaptation mechanism to capture the achievable rate. It is 
made possible, since our assumption on bounded noise rise 
density allows the computation of the uplink transmission 
power density regardless of the schedule. The transmission 
power density is sufficient for the rate adaptation. Clearly, in 
practice, a dynamic rate adaptation mechanism would better 
estimate the achievable rate and facilitate a better utilization 
of the multi-user diversity resulting in a higher throughput. 

VII. Simulation results 
A. Numerical Results - Shannon Capacity 

In this sub-section, we give the numerical results for the 
algorithms in section [V] and section I VII and compare them to 
a fixed power scheme. Our main figure of merit is the total 
throughput in the system, summed over all frames and all 
cells. However, to stress out the benefits of the noise-rise based 
schemes, the standard deviation of the interference seen at 
the cells as well as the actual power used are also given. In 
particular, the following simulation was carried out. 

The considered deployment include 72 hexagon-shaped 
cells (base stations), with 722 users (mobile stations). The 
user locations were drawn uniformly at random on the space, 
with the restriction of a minimum of 2 users per cell. The 
topology thus created is given in Figure Q] The path losses (in 
dB) were calculated according to the COST-Hata model [19 



Capter 4, equation 4.4.3] 

PL =46.3 + 33.91og 10 (/) - 13.82 log 10 (M - ah m 

+ (44.9 - 6.55 \og 10 {h b )) log 10 (d) + c m , 

where / = 2000MHz was the transmission frequency, hs = 
50 was the base station antenna effective height, ah m was the 
mobile station antenna height correction factor and c m was 
taken as for a medium size cite. The simulation included 80 
frames, with the topology and path gains fixed throughout. 

The algorithm used to solve the optimal, noise-rise con- 
straint problem (fT5l i was Algorithm Iterative-Water- 
Filling. In the noise-rise density constraint algorithm only 
a single user per cell is scheduled according to Algorithm 
Resource-Allocation. That is, this user is selected ac- 
cording to the maximal value of u>iXi log ^1 + , assuming 
the user exhausts the noise-rise budget on its own. In the 
fixed power scheme, again, only a single user per cell is 
scheduled. However, in this schemes, the allocated power 
is fixed at the same constant P for all users, and only the 
user which under such an allocation will have the maximal 
normalized rate cjj log fl + "^r 1 ) is scheduled. Of course, for 
a fair comparison, the fixed power P is set such that the mean 
interference caused by all schemes is the same. 

The weights cjj were initialized to a constant value, and 
were updated after each frame according to ifTTI 

T i {t)=T i (t-l) + (l-/3)-B i {t-l) 

and 

Ui (t) = 1/T<(t), 

with (3 = 0.9 and Bi(t — 1) is the number of bits delivered 
to user i at sub-frame £ — 1 (0 if user i was not scheduled at 
frame t — 1). 

In all three schemes, the actual throughputs after each frame 
were calculated according to the true SINR at the receivers 
(Shannon capacity for band-limited Gaussian channel). This 
is to simulate the exact scenario for which out analytic claims 
apply. Note, however, that this is fundamentally different from 
the all-encompassing simulation described in the next section, 
where specific modulations and packet loss rates are taken into 
account. 

The results are given in Figures [2] to [4] In Figure |2j the total 
throughput, over all cells and frames is given. The benefit of 
both noise-rise schemes over the fixed power scheme is clear 
for all noise-rise values (equivalently, all SNR values). The 
benefit of the optimal algorithm over the sub-optimal is also 
clear, as it indeed usually schedules more than one user per 
cell. Note that it is possible, though, to construct topologies 
in which the benefit is small. 

A key advantage in the noise-rise schemes is that the 
interference seen at the neighboring cells is concentrated 
around the fixed noise-rise value /, allowing them to plan 
modulation and coding accordingly. Clearly, when the variance 
of the actual interference observed at the receivers is high, it 
is harder to choose the appropriate modulation and coding, 
forcing the senders to either aim at lower modulation and 
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Fig. 1. The network topology. 72 cells with 722 mobile stations, distributed 
uniformly at random. Each cell was restricted to have a minimum of 2 mobile 
stations. 



coding schemes or suffer a high packet loss. The variance of 
the observed interference for all schemes is given in Figure [3] 
Note that the average is fixed to the same value in all schemes. 

Finally, we address the lack of maximum power constraint 
(head room). The four figures in Figure|4]depict the histograms 
of the actual powers allocated to the users in the optimal noise- 
rise scheme. It is clear that a substantial fraction of the users is 
allocated powers below or around that allocated by the fixed 
power scheme, and even the small fraction with high powers 
does not require a consequential increase in power. Thus, 
although the optimal scheme is not head room constrained, 
its actual power usage is moderate. 

B. An IMT-Advanced Simulation 

We illustrate the advantages of the noise rise approach 
via extensive system simulations. To that end, we developed 
a Matlab simulation that comply with the guidelines for 
evaluation of radio interface technologies for IMT-Advanced 
0. Specifically, we simulate a two tier hexagonal deployment 
with 19 sites, each contain three sectors. Each sector operates 
as an independent base station. Users are located uniformly 
across the deployment with 10 users in each cell. Statistics 
are collected from all cells, which is facilitated by a wrap 
around geometry such that all cells suffer from interference 
from two tiers neighboring cells. 

Following the IMT-Advanced guidelines, we implemented 
the WINNER II |20| stochastic channel modeling. The model- 
ing approach is based on the geometry of the network layout. 
The large-scale parameters such as path loss and shadowing 
factor are generated according to the geometric positions of the 
base station and the user. Then the statistical channel behavior 
is defined by distribution functions of delay and angle as well 
as by the power delay and angular profiles. We consider the 
urban macro-cell scenario, where typically mobile stations are 




Noise Rise in dB 

Fig. 2. The commulative throuhput over all frames (80) and all cells. The 
optimal algorithm of section|V]is compared to the sub-optimal algorithm given 
in section IVTl and the fixed power scheme. 
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Fig. 3. Standard deviation of the interference seen at the cells (over all 
frames and all cells). The optimal algorithm of section [V] is compared to the 
sub-optimal algorithm given in section [VTl and the fixed power scheme. While 
the average interference is constrained to the same value in all schemes, the 
two noise-rise schemes (the optimal and the sub-optimal algorithms) exhibit 
a significantly lower standard deviation. That is, the interference seen under 
these schemes is centered around the average value, allowing better rate 
planning by the base stations. 
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Fig. 4. Histograms of the allocated powers in the optimal algorithm (only 
powers above 0.1 are included). The dashed red line represents the power 
allocated by the fixed power scheme. It is clear that the number of users 
with high powers compared to the fixed power scheme is small. Moreover, 
even those users do not use more than double or three times as much power 
compared to the fixed power scheme. 



located outdoors at street level and base stations are fixed 
clearly above surrounding building heights. The assumed inter- 
site distance is 500m. We assume two receive antennas at the 
base station and one transmit antenna and the mobile station. 
The maximal uplink transmit power is 24dBm. The user and 
the base station antenna gain is and 17 dBi, respectively. The 
assumed noise figure at the user terminal and the base station 
is 7 and 5dB, respectively. 

We consider the IEEE 802.16m uplink frame structure 
(Time Domain Duplexing - TDD mode). The assume total 
bandwidth is of 10MHz, which is occupied by 48 resource 
units. The simulation runs sub-frame by sub-frame, performing 
the scheduling and power allocation at the beginning of each 
sub-frame. We run the simulation over 1000 sub-frames for 
5 drops. In each drop all users are reallocated across the 
deployment. 

The exponential-effective SINR Mapping (EESM) approach 
ll2~T1l is used to map the system level SINR onto the link level 
cursves to determine the resulting block error rate. Upon block 
errors, a synchronized retransmission is scheduled with the 
highest priority after 5 sub-frames from the original transmis- 
sion. Upon reception of retransmission the receiver performs 
chase combining. The maximal number of retransmissions is 
4. 

We assume a full-queue model where all users have back- 
logged traffic. The weights uii were initialized to a constant 
value, and were updated after each frame according to Ti (i) = 

Ti(t-l) + (l-P)-Bi(t-l) and^(i) = 1/T t (t), with various 
values of /3, 

Figure [5] compares three power control schemes, namely, 
(i) constrained noise rise (termed N.R) based on Algorithm 
Iterative-Water-Filling; (ii) constrained noise rise den- 
sity (termed N.R. density) based on Algorithm RESOURCE- 
ALLOCATION; and (iii) target received SINR (termed SINR). 
The last scheme represents the traditional power allocation 
scheme (which aims at obtaining the required SINR level 
at the receiver) and is used as a benchmark for our noise 
rise approach. One can see that both noise rise algorithms 
obtain a dramatically better spectral efficiency. For example, 
the constrained noise rise approach obtains spectral efficiency 
of 1 [bit/Hz/Sec] for (3 = 0.99 where the SINR based scheme 
obtains only 0.5 [bit/Hz/Sec]. Yet, the traditional approach 
obtains a slightly higher cell edge spectral efficiency. Inter- 
estingly, the constrained noise rise density algorithm obtains 
a better spectral efficiency than the original scheme. This 
advantage is due to the difference between the Gaussian Mul- 
tiple Access Channel model and the simulated IMT advanced 
EESM channel model, which does not disrupt the constrained 
noise rise density approach (as discussed in section [VTb . 

VIII. Conclusion 

In this paper, we considered a joint scheduling and power 
allocation problem. Specifically, to mitigate inter-cell inter- 
ference, we suggested a novel approach, which considers the 
interference caused by a base station to neighboring base 
stations as a resource to be allocated, similar to bandwidth 
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Fig. 5. Uplink spectral efficiency (a) average (b) cell edge. 



or power. The essence of this approach is as follows. Each 
user, based on its channel gains, creates a different interference 
to neighboring cells when transmitting. A base station, while 
allocating power and bandwidth to its subscribers, does so 
in a way such that the total aggregated interference its users 
generate does not pass a certain value - its noise rise budget. 

We rigorously formulated the problem as a convex opti- 
mization problem with linear constraints and suggested an 
efficient iterative algorithm for its solution, based on known 
and new water-filling based solutions to its separate problems. 
We then devised a modified algorithm, which optimizes power 
and bandwidth allocations assuming the noise rise is constraint 
for each sub-channel (a portion of the bandwidth allocated to 
a single user). This algorithm is highly efficient, yet performs 
better than known techniques, approximating the performance 
of the optimal algorithm. These performance guarantees where 
also depicted in two extensive simulations, one based on the 
analytical expressions and one based on practical system built 
based on the IMT-Advanced specification. 

Potential applications of the suggested algorithms include, 
but are not limited to, uplink scheduling in IEEE 802.16e/m, 
3GPP LTE, and more. 
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Appendix A 
Proof of Proposition[T] 

Proof: Consider the negative utility function 

- ^2 UiXi log ( 1 + J . (26) 

We first show that this function is convex. To this end, consider 
a single summand 



f(xi,pi) = -uj l x t log ( 1 + 



The Hessian matrix, restricted to the variables x.- L and pi is 
given by 

/ df _df_ 





Note that 




ai- 



Oi 2 



Hence, the Hessian matrix of (|26l , H(x,p), satisfies 
(ai, . . . ,a 2 M)H(x,p)(a-L, . . . , a 2 Ai) T 



E 



(1 



otii-i— - u 2 i I > (27) 



II 



for all (ai,...,0!2Ar) G 



p2M 



Appendix B 
Proof of Proposition!!] 

Proof: In Appendix [A] we show the convexity of the 
optimization problem. The proposition will now follow from 
a straight forward application of the KKT conditions ifTTl 
Section 5.5.3]. Namely, we write 

^2 U i X i lo S ( 1 + ) + ^1 liVl ~ 1 

teM \ x< ' \ieM J 

+ 1 Xj - 1 j - ^iXi - (28) 

and the proposition follows by requiring VL = and that for 
all i G A4 we have 

/iift = 0, fiiXt =0, /ij > 0, \Xi > 0. 



Appendix D 
Proof of Corollary Q] 

Proof: We first show that the negative utility func- 
tion — YlizM 1°§ y- j satisfied an existence and 
uniqueness constraint [18 Section 2]. That is, fixing any 
2M — 1 variables and optimizing on the remaining variable, 
the resulting problem has a unique (global) minimizer in the 
range. This is easily seen from the Hessian matrix calculated 
in the proof of Proposition as the negative utility function 
is convex in each of the variables. 

We now mention that the transform (p, x) t = T[(p, x) t -i] 
defined by one iteration of Algorithm ITERATIVE- Water- 
FILLING has no fixed points (for which T[v] = v) besides 
the global optimum of ( TToT ). This is, again, since the negative 
utility function is convex. The corollary will now follow by 
applying [18, Theorem 2]. ■ 



Appendix C 
Proof of Proposition[3] 

Proof: To prove item 1) we proceed as follows. Set = 
a and consider the function /(a) = log(l + a) — It easy 
to verify that /(0) = and that f(a) is non-negative and 
monotonically increasing for any a > 0. Thus, for any ^ < 0, 
the equation f(a) - 



V i e 



— will have a unique solution at some 
- > 0. Hence, Xi > 0. Moreover, for fixed pi,ei,u!i, 
the smaller A2 is, the larger is the solution to f(a) = — 
that is, the smaller ii, for all i, and hence the smaller is J^. Xi. 

To prove item 2) we note the following. Since . ii is 
monotonically increasing in A2, the value of A2 such that 
= 1 i s c l ear ly upper bounded by the value of A2 for 
which the "weakest" user, the user which requires the largest 
A2 in order to achieve Xi = 1, indeed gets it (since in this 
case all other users will have x i< > 1 and ^ . ii will clearly 
surpass 1). Thus, the optimal A2 is at most 



max {A2 



s.t. 



= 1} = 



max u)i 



1 + Piei 



log(l +Piei) 



(29) 



On the other hand, the value of A2 such that J^i %i = 1 
is clearly lower bounded by the value of A2 for which the 
"strongest" user, the user which requires the minimal A2 in 
order to have Xi — 1/M, indeed gets it (since in this case 
all other users will have ii> < 1/M and thus Y^i&i W 1H ^ e 
strictly smaller than 1). As a result, the optimal A2 is at least 



min {A2 

ieM 



s.t. 



= 1/M} 

1 + M Pl ei 
This completes the proof. 



mm < 0J4 



log (1 + M Pi ei) 



(30) 



